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ABSTRACT 


A software package was developed to provide the inputs for a 
passive single-sensor acoustic tracker. This required a computation 
of sound propagation velocity and resolving high-resolution frequency 
information from raw acoustic data. A selected frequency was then 
automatically tracked from an up-doppler to a down-doppler condition. 
This frequency information coupled with associated bearing information 


was then used to provide continuous target fixes. 
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I. INTRODUCTION 


Continuous fixing of an underwater sound source by passive means 
has traditionally been carried out by using multiple sensors. Single 
sensor fixing requires varying power coefficients, frequencies, and/or 
bearing changes and usually produces an estimated position with an 
accuracy that varies greatly with the type of sensor and processing 
utilized and the skill of the evaluator. Anaccurate single sensor track 
would be advantageous in that positional errors could be minimized with 
only one buoy or array required, providing greater flexibility in that a 
series of sensors are not required to be monitored, and a decrease in 
buoy expenditures could be expected. To provide the necessary frequency 
data, 2 software package was developed that computed the sound propaga- 
tion velocity, detected a discrete frequency signal and automatically 
tracked that signal from a full up-doppler through a full down-doppler 
situation. This frequency output combined with associated bearing in- 
formation, and a value for the sound propagation speed with an estimate 
of the speed of the target provides sufficient input data for single sensor 


fixing. 
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Il. THE TRACKING PROBLEM 


A. STATEMENT OF THE PROBLEM 

Chapter 3 of [1] developes the required input quantities for passive, 
single-sensor fixing. These values are an estimate of target speed, the 
velocity of sound propagation in the medium, the doppler shifted fre- 
quency and a bearing from the sensor to the target. The estimate of 
target speed would have to be determined by tactical considerations at 
the time and location of signal acquisition, i.e., target transiting or on- 
station, signal strength, detection ranges, etc. ‘The bearing would be 
provided by a standard DIFAR buoy or anarray. In this problem the 
standard deviation of bearing measurement noise was varied from Bo Lo 


oO 


30°. Methods of providing the sound propagation velocity and the doppler 


shifted frequency were then devised. 


B. DEVELOPING THE INPUTS 

Equation 1.1. of [2] provides a good approximation of sound speed 
in sea water as a function of temperature, salinity and pressure. It is 
repeated here as equation l. 

c(T, S, z)=1449+4, 62T-. o54T 41. 3(S-35)+. 0172 (1) 
where c is in meters per second, T is in degrees centigrade, S is in 
parts per thousand, and zis in meters. Converting this equation so 
that c is in yards per second, T is in degrees Farenheit and depth is in 


feet, yields equation 2, 
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If the target depth is known, the median depth of the target and hydro- 
phone should be used. Otherwise, use the hydrophone depth. Salinity 
would normally be set to 35ppt unless operating in water of a known 
different salinity. The resultant output from this equation would be set 
equal to the velocity of sound propagation, VP, and used as an input 
value to the target tracker program. 

The frequency requirements called for a high-resolution value 
which dictated a relatively long time window. A twenty second window 
was used which resulted ina .05Hz FFT resolution. This resolution 
was further enhanced by interpolation using the relative values of the 
power coefficients associated with the discrete frequency being tracked. 
In order to maintain contact on the same discrete frequency line, a 
frequency tracker algorithm was devised. This provided the necessary 


frequency data for the target tracker program. 
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Ill SOFTWARE IMPLEMENTATION 


A. SOUND SPEED 

This process was a simple matter of implernenting the equation for. 
sound speed in seawater, equation 2. Providing the chosen depth, ob- 
served temperature, and salinity results ina best estimate of the actual 
sound speed in the operating area. Updating this value would be necessary 


as the operating area or environmental conditions changed. 


B. FREQUENCY 
1. Parameters for Frequency Computations 

In order to provide realistic frequency inputs, a magnetic tape 
was obtained that contained analog, raw acoustic signals representing 
arnbient noise plus an underwater sound emitter with various discrete 
frequencies. The doppler shift of a discrete frequency signal is directly 
proportional to the rest frequency of the signal. Thus doppler shift 
measurements are more easily obtained from higher frequencies. 
However, signal attenuation in the medium is much higher at higher 
frequencies. At present, frequencies in the 150Hz to 300Hz range will 
yield usable values with realizable frequency resolution measurements. 
The highest signal frequency of interest determines the required sam- 
pling rate of the data. The sampling rate and the time window length 


determine the required dimension of the data input matrix. With these 
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considerations in mind, the following parameters were chosen in digitiz- 
ing the data: 
‘Sampling Rate of 512 samples/second yields Nyquist frequency 
of 256Hz; 
Input data low-pass filtered to 250Hz to prevent aliasing errors; 
Time record length of 20 seconds for FFT resolution of .05 
Hz/Bin; 
Hanning window function applied in time domain to reduce 
sidelobes associated with rectangular data windows. 
2. Power Spectrum 
With the above chosen data parameters, the input matrix to the 
POWER Subroutine™ [3] was dimensioned to 20,480 bins. This cor- 
responds to two bins per data point (1 cosine termand 1 sine term for 
10,240 data points) or 10,240 Complex Fourier coefficients. The real 
data values were loaded into the odd (cosine or real) bins and zero's 
were provided for the even (sine or imaginary) bins. On return from 
the FFT Subroutine [6] to the POWER Subroutine, the cosine term is 
squared and added to its associated sine term squared thus representing 
the power at some discrete frequency. This discrete frequency is an 
integer multiple of the FFT processing resolution (the reciproca] of the 
time record length). The first bin, Y(1), contains the zero frequency 
power coefficient, the second bin contains the .05 Hz power coefficient, 
and so on in .05 Hz steps. Thus bins Y(1) through Y(5001) now contained 
the power coefficients from zero to 250 Hz. 


“Subroutine STATS and POWER were written by Arfman, J. F., Jr., 
and were adapted for use in this program. 
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3. Search Window and Background Noise Suppression 

With over five thousand power coefficients available for track- 
ing, it is convenient to limit the range of search to those that offer the 
best data for doppler tracking. A signal having favorable characteristics 
is first identified, i.e., frequency above 100 to 150 Hz, apparently 
stable on a one Hz resolution processor, and of sufficient strength to 
maintain a favorable aigaaileglatiee ratio for a stash of the time. 
Then a search window is designated by the adjustable parameters FL 
and SRBW, where FL is the lower frequency limit of the search window 
and SRBW is the search window bandwidth. Both parameters must be 
some integer multiple of the FFT resolution and would usualiy be kept 
at whole numbers for convenience. The total number of power coefficients 
in the search window would equal N=20(SRBW) and are designated YA(1) to 
YA(N). 

Ideally when the frequency to be tracked (FTRK) has been iden- 
tified, FL and SRBW should be chosen such that FTRK is equalto FL 
plus (1/2)SRBW and SRBW is as large as is reasonable to keep the ratio 
of signal power in the window to the total power in the window much iess 

L M 
than one. If rae YA(K) KK yi YA(J) for L signal bins and M noise 
K=1 J=1 
bins with L+M=N, then the summation as I goes from one to N of YA(I)/N 
is a reasonable approximation of the ambient noise level around FTRK [4], 
[5]. This calculation is performed and set equal to YNOIS. Then YNOIS, 
or some function of it dependent on the signal-to-noise ratio, is sub- 


tracted from each power coefficient in the search window. Under the 
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assumption that the noise power spectrum is constant (white noise) 
within SRBW, only the signal power spectrum remains in the window. 

4, Frequency Tracking to Follow the Doppler Shift 

The frequency selected to be tracked will be the strongest sig- 

nal in the search window when the first transform is computed. The 
bin number of that signalis saved as TBIN(JJ). In each succeeding 
transform the strongest signal (largest coefficient) in the search window 
is designated TBIN(JJ) and is compared to TBIN(JJ-1). A bin shift 
greater than 10 is determined to be a new or different discrete signal. 
The new power coefficients around the last center bin are then examined 
for discrete signal energy. The maximum value found would be con- 
sidered the signal energy being tracked if its bin number is equal to the 
last center bin number plus or minus Il, and its magnitude is greater 
than a chosen threshold minimum. In this problem this value was set 
to three times the ambient noise value. If the maximum coefficient 
found is more than one bin from the last center bin, the old center bin is 
examined. If its coefficient is larger than at least one of its sidelobe 
bin values and at least half as large as the maximum found and larger 
than the chosen threshold minimum value, the old signal bin is con- 
sidered to still contain the signal being tracked. If any of these test 
fail, the new maximum found is cae see to the chosen minimum 
value. If the new maximum is larger, it is accepted as the signal being 
tracked; if not the contact is considered lost. Figure 1] possibly 


clarifies the above process. 
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LOAD NEW POWER COEFFICIENTS 
OF BINS AROUND TBIN(JJ-1) AND 

| FIND NEW MAXIMUM COEFFICIENT. 
RESET TBIN(JJ) TO NEW VALUE AND 
COMPUTE NEW BNSHF. 


>) BNSHF > 10 | 


-. z — on — . F 
> BNSHF S 1 PP MAX COEFF> THMIN fF 


At 


| MAX COEFF IS 
SIGNAL BEING 
TRACKED 


[CHECK NEW COEFF 


Sie pBNSHr=10 °° Lo) 
is ao OF TBIN(JJ-1) 


; F 
ae ae Ee ARTO 


amon dey 

( COEFF(TBIN(JJ-1))> COEFF((TBIN(JJ-1))+1) Jooend 
PCR AL t: es i BY 

[ COEFF(TBIN(JJ- i) > (COEFF(TBIN(JJ)))- 1/2 } 


( COEFF(TBIN(JJ-1))> THMIN } tT 
DTT EIE TA SF 


arte. 


A ARTE, 


OLD TBIN CONTAINS 
SIGNAL TO TRACK 


COEFF(TBIN(JJ))> THMIN 


BIN SHIFT = BNSHF = TBIN(JJ) » TBIN(JJ-1) 
THMIN = A CHOSEN THRESHOLD MINIMUM VALUE 


FIGURE 1. FREQUENCY TRACKER DECISION FLOW. 
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Once the bin containing the largest portion of the signal energy 
being tracked has been established, the frequency of that signal is 
refined by computing the first moment of the three bin power herent 
around it. If the detected signal is an integer ae of the EFT 
resolution, adjacent bin spill over will be minimum and the refining 
calculation would make no correction to the frequency of the bin 
associated with the maximum power. When the detected signal is not 
an integer multiple of the FFT resolution the refining calculation adjusts 
the frequency to a mid-resolution value corresponding to the first mo- 
ment of the hump. See Figure 2. The frequency corresponding to bin 
number K is stored inthe YL matrix as YL(K). The frequency refining 
equation used was 


FTRK(JJ)=YL(K)+RESOL*(YA(K+1)-YA(K-1))/(YA(K41) 
+YA(K)+YA(K-1)) (3) 


where 
FTRK=the discrete frequency being tracked, 
RESOL=the resolution of the FFT processing, 
YA(M)=the power coefficient of the Mth bin. 
Letting YL(K)=220.25Hz and using the sample values from figure 2, 
FTRK(JJ) = 220.250Hz + .05Hz - ( as seen ans 
= ZL0i250Hast .05Hz ~ (=. 25.) 
= 220.250Hz - -0125T1z 
= 220.2375Hz 


220.238Hz (rounded) 
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@--— FTRK=220, 238Hz 
“\ 
I teGem YL (K)=220. 250Hz 


YA(K) 
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RESOLVING MID-RESOLUTION FREQUENCIES WHEN THE 
DETECTED SIGNAL IS NOT AN INTEGER MULTIPLE OF 


FIGURE 2. 
THE FFT RESOLUTION. 
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Without this calculation the returned value of the frequency being tracked 
would have been 220.250Hz. A five bin computation was found to offer 
no significant improvement over the three bin computation due to the 


increased probability of noise or adjacent discrete signal influence. 


C. DESIGNATING PROBLEM VALUES 
In computing the sound propagation speed for the following examples, 
nominal values were used because the true values were not available. 
Values used were: 
Temperature = 48 F 


Depth =O UREL 


Salinity 35 ppt 

The sound emitter on the data tape appeared to make a slow turn 
into the buoy for the first 4 minutes of the tape and then maintain a 
constant course at a speed of approximately 4.5 knots. The closest- 
point-of-approach (CPA) to the hydrophone was approximately 300 to 700 
yards and 20 minutes of data was utilized with CPA at about the 11 to 12 
minute point on the tape. A 20 Hz search window bandwidth was used 
with several different lower frequency limits on the window. It was 
found that one transform per minute comfortably contained the tracked 


frequencies within the 21 bin tracking window in all cases for a slow- 


moving, near-field emitter. 
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Pee DATA OUPRPUT 
1. Run 1 (See Table 1) 
With FL set to 40 Hz and SRBW at 20 Hz, the rest frequency 
of the tracked signal was approximately 50.13 Hz. The total doppler 
shift was only 0.15 Hz and was too low to provide good tracking data. 
2. Run 2 (See Table 2) 
For this run FL was set to 220 Hz with SRBW still at 20 Hz. 
Two strong discrete frequency signals were within the search window 
bandwidth, one's rest frequency at 233.87 Hz and the other at 223, 22Hz. 
The maximum power coefficient changed between these two signals 
several times. In the initial transform the tracker locked on to the 
lower of the two frequencies and a continuous track was maintained even 
though the signal faded out and then came back. This run gave the fre. 
guency tracker a good check-out in that (a) the maximum coefficient in 
the search window had considerable variations between discrete signals 
(b) a noise spike occurred in the tracker window at time=430 (FILE=220) 
and (c) the tracked signal faded completely out at time=1090 (FILE=550) 
and was regained one minute later. The total doppler shift was approx- 
imately 0.6 Hz and provided good data for doppler tracking. The computer 
output for this run is included in this report in its entirety. 
3. Run 3 (See Table 3) 
For this run FL was set to225 Hz and SRBW remained at 20 Hz. 


The only significant signal in the search window was the 233. 87 Hz signal 
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Time Center Ratio: Frequency 


(seconds) Bin Signal/Noise Tracked 
10.0 204 155/66 50.155 
70.0 205 1759/38 50.190 

130.0 205 1619/59 50.185 
190.0 205 2001/54 50.203 
250.0 205 1620/39 50. 180 
310.0 205 1527/49 50.192 
370.0 204 878/40 50.197 
430.0 205 247/42 50.151 
490.0 205 9320/78 50.178 
550.0 205 13280/68 50.191 
610.0 204 56746/254 50. 168 
670.0 203 38878/144 50.105 
730.0 203 42170/187 50.085 
790.0 203 47640 /227 50.081 
850.0 202 6941/55 50.060 
910.0 202 9654/71 50.063 
970.0 202 18841/125 56.053 
1030.0 202 17558/111 50.057 
1090.0 202 17546/104 50.054 
1150.0 202 14377/87 50.050 


Rest Frequency approximately 50.13 Hz 


TABLE I. 
FREQUENCY TRACKER OUTPUT WITH FL= 50 Hz, SRBW=20 Hz. 
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Lime Center Ratio: Frequency 


(seconds) Bin Signal/Noise Tracked _ 
10.0 70 69/6 223.437 
70.0 71 112/5 223.481 

130.0 71 61/5 223.495 
190.0 72 52/5 223.556 
250.0 71 254/6 223.489 
310.0 70 218/6 223.473 
370.0 69 490/9 223.420 
430.0 69 64/8 223.407 
490.0 70 743/24 223.434 
550.0 71 2131/41 223.488 
610.0 69 10429/68 223.396 
670.0 62 4795/50 223.063 
730.0 60 2554/26 222.960 
790.0 61 2119/16 222.995 
850.0 60 322/14 222.938 
910.0 59 245/10 222.882 
970.0 59 60/4 222, 898 

1030.0 59 112/6 222. 886 

1090.0 60 12/5 Lost Track 

1150.0 59 115/5 222.901 


Rest Frequency approximately 223.22 Hz 


TABLE II. 
FREQUENCY TRACKER OUTPUT WITH FL=220 Hz, SRBW=20 Hz. 
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Time Center Ratio: Frequency 


(seconds) Bin Signal/Noise Tracked 
10.0 183 61/4 234.100 
70.0 184 42/4 234.151 

130.0 185 21/4 234.193 
190.0 185 vara Ae) 234.224 
250.0 184 167/5 234. 169 
310.0 184 160/5 234. 163 
370.0 183 829/7 234.102 
430.0 184 708/8 234.144 
490.0 184 3295/19 234.132 
550.0 184 4824/32 234. 130 
610.0 183 4669/30 234.079 
670.0 176 6049/29 233.756 
730.0 174 906/13 233.644 
790.0 174 766/7 233.634 
850.0 172 1082/11 23375105 
910.0 172 1020/8 233.558 
970.0 172 27/4 235-5103 

1030.0 172 357/5 233.545 

1090.0 172 135/5 23 3e 5310 

1150.0 7A 23/5 233.525 


Rest Frequency approximately 233.87 Hz 


TABLE III. 
FREQUENCY TRACKER OUTPUT WITH FL=225 Hz, SRBW=20 Hz. 
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and this was tracked with only a minimal amount of decision require- 
ments around CPA. The total doppler shift was approximately 0.65 Hz 


and provided good data for doppler tracking. 


E. TARGET TRACKING 

Since no bearing information was available corresponding to the 
frequency data, the IBM 360 Computer was used to generate simulated 
real time bearings. This was accomplished by programming a pseudo 
"true track" into the computer so that the target speed and the CPA time 
were matched as wellas possible. The initial starting position was 
X=500 yards, Y=-1050 yards at Ene T(1)=190 seconds. The heading 
and speed used were 90 degrees and 4.5 knots where zero degrees 
corresponds to the positive X axis. The computer then computed "true 
bearings'' to the target and added random noise values of a programmed 
standard deviation, SA, to these values. Thus the true track plotted 
by the computer is only a best estimate of the targets relative position 
from the buoy. 

The X-Y Filter developed by Mitschang [1] was used to process the 
data and develop the target track. To initialize the filter and commence 
tracking, the measured bearing must have changed by at least three 
standard deviations of the measurement noise and the doppler shifted 
frequency must have decreased from its value at the start of the initializa- 
tion process. A modification to the filter was attempted so that while 


initializing, an increasing frequency would be detected as a maneuvering 
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target. Time, bearing, and frequency values occurring with the maxi- 
mum frequency detected would then be retained as the initial values to 
be used when the bearing shift and non-increasing frequency criteria 
were again satisfied. This would delay the start of the tracking process 
but would prevent data from a maneuvering target being used to initialize 
the state equations based ona constant heading, constant speed target. 
However, indexing problems were encountered and ae not resolved 
in time to include the modification in this report. 

The following difference equations developed in Chapter 3 of [1] 


were used to compute the initial range, heading, and rest frequency of 


the target. 
Re e(VEADRne DELT (asta (DELeI" © 5 (4) 
@si=Qave-180 - arcsin((-VP: DELF)/(f£(last): DEL@- Vsi)) (5) 
F0i=Fave [1 - (vsi/VP)cos(@si - 8ave)] (6) 
where 
VP =speed of sound propagation 


DELF =doppler shift of the frequency 


DELT =time required to meet initialization criteria 
{f(last) =last frequency measured for initialization 
DEL®# =change in bearing from sensor to target 

6si =initial estimate of target heading 

Bave =average bearing for the initialization period 
Vsi =initial estimate of target speed 
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Foi =initial estimate of the rest frequency 
Fave =average frequency for the initialization period 
Then based ona constant heading, constant speed target, the initial 


filter state equations are: 


X(1) =X position = R cos Save (7) 
X(2) =X velocity = Vsi cos Osi (8) 
X(3) =Y position = R sin Oave (9) 
X(4) =Y velocity = Vsi sin Osi (10) 
X(5) =Rest frequency = FOi Ci) 


The initial covariance matrix values are computed by the "direct 
method'' developed in Chapter 4 of [1] and an extended Kalman filtering 
technique was used to refine the estimates of the target's position, 
heading, speed, and rest frequency as the problem progresses. 

Since the target was maneuvering for the first four minutes of the 
tape, initialization of the target tracker was commenced at time 
T(4)=190 seconds, with the values from Table 3 as input data. The 
average initialization time, AVKJ, varied with the standard deviation of 
the bearing measurement noise, SA. This time varied from about 4 


minutes for SA=5 degrees to about 8 minutes for SA=30 degrees. 
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IV. RESULTS OBTAINED FROM REAL FREQUENCY INPUTS 


Once the frequency tracker has locked on to a stable but doppler 
shifted frequency, the inputs required of the operator are an estimate of 
the rest frequency, FO, and an estimate of the target's speed, AVSK. 
‘These are "priming" values used in computing the initial heading of the 
target and affect the ultimate accuracy of the computed track. A third 
variable affecting the initialization is the standard deviation of the bear-~- 
ing error, SA. As SA increases, the filter takes longer to start tracking 
and has more data on which to base its original estimates. 

The program was run for various combinations of the three variables 
FO, AVSK, andSA. The true vaiue for FO was approximately 233.87 Ma 
and for AVSK was 4.5 knots. An estimated FO greater than the detected 
frequency causes a "down-doppler" solution for the initial course. Con- 
versely an estimate of FO less than the detected frequency results in an 
"up-doppler" solution for the initial course. Each combination of var- 
iables was looped so that one hundred runs were made with that set of 
inputs and the averaged filtered positions were plotted. The ''+'' symbol 
on each plot represents the computed position of the target at time T(1) 
plus one-half of the initialization period. The range is computed by 
Equation (4) and the bearing is the average bearing for the initialization 


period. The true points from which the noisy bearing measurements 
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were made are plotted with the symbol '*'"'. Averaged filter positions 
are plotted with an 'X'". The ascending order of precedence on each 
plot is +,*,X. Figures 3 through 24 show the performance of the system 
with real, high resolution frequencies as inputs. F(DET) signifies 


initial frequency detected. 
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FIGURE 3. FILTER OUTPUT WITH TRUE FO AND AVSK INPUTS, 
ZERO BEARING ERROR, AND TEN MINUTE INITIALI- 
ZATION PERIOD. FIFTH TRUE POINT AND FINAL 
TRUE POINT SUPPRESSED BY FILTER OUTPUT 
POSITION 
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HIiGURE 4, FILTER OUTPUT FOR EST. FOS FH(DET), EST. 
AVSK * TRUE AVSK, AND SA = 30°. FINAL TRUE 
POINT SUPPRESSED BY FILTER OUTPUT FINAL 
POSITION. 
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FIGURE 5. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER THAN 
F(DET) EST. AVSK < TRUE AVSK, and SA = 30°. 12th 
and 13th TRUE POINTS SUPPRESSED BY FILTER OUTPUTS 
13th and 14th POINTS. 
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FIGURE 6. FILTER OUTPUT FOR EST. FOS F(DET), EST. AVSK> 
TRUE AVSK, AND SA = 30°. 


SS 


wer 


sby 0001 - 


1000 yds - 


vk 


r 


mes career a a eS 


-1000 yds - ¥ 


0 1000 yds 


FIGURE 7. FILTER OUTPUT FOR EST. FO> F (DET) EST. AVSK> 
TRUE AVSK SA = 30°. 


34 


Aoa TUSTUO ARTS 
‘or « AG SEVA aU3BT 


FIGURE 8. 


1000 yds 


-1000 yds 


BRILTER OUTPUT FOR EST: FOS 


aaa a fe 


1000 yds 


Dp) 


TRUE AVSK BY 90 PERCENT, SA = a0)e, 


3D 


> 


EST. AVSK> 


abe O00! . 


ity vi 0 

: oo ACL TUNTUO AAT aS 
(hi THEO 00 ¥4 ARBVARUAT 
6 
# 
, i 


1000 yds + 


i 
i 
apc apt e Siearet res aon pss oer eon rarer th arn eraeer nap esa a oer oes 
! 
| 
} 
| 
} 
or 
bq 1 
| 
Pee een Herren ee Bet Ue Oh eve eee 


Bes 


ie 


PNOOU Cs bora ee 


FIGURE 9. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F(DET), EST. AVSK< TRUE AVSK, AND 
SA = 15°. 12th and 13th TRUE POINTS SUPPRESSED 
BY 13th and 14th FILTER OUTPUT POINTS. 
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FIGURE 10. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F(DET), EST. AVSK % TRUE AVSK, AND SA = 
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FIGURE 11. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F (DET),EST. AVSK > TRUE AVSK, AND SA = 
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FIGURE 12. FILTER OUTPUT FOR EST. FO > F(DET), EST. 
AVSK > TRUE AVSK, SA = 15°. 
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AVSK > TRUE AVSK, SA = 15°, 
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FIGURE 14. FILTER OUTPUT FOR EST. FO=> F(DET), EST. 


AVSK > TRUE AVSK, SA = 15°. 


4] 


: Se eee 
‘oR .Tea ACT TOITUO ABTA 
) 2 HEVA YUASA <MRYA 


(ooGhyds =) 


NOOO) yds te 


0 1000 yds 


FIGURE 15. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F(DET), EST. AVSK< TRUE AVSK, AND 
SA = 10°. 11th TRUE POINT SUPPRESSED BY 11th 


FILTER OUTPUT POSITION. 
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FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F(DET), EST. AVSK&¢ TRUE AVSK, AND 
SA = 10°. 12th THROUGH 16th TRUE POINTS 
SUPPRESSED BY 13th THROUGH 17th FILTER OUT- 
PUT POSITIONS. 
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FIGURE 17. FILTER OUTPUT FOR EST. FO< F(DET), 
EST. AVSK > TRUE AVSK, SA = LO 
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FIGURE 18. 
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FILTER OUTPUT FOR EST. FO> F(DET), 
EST. AVSK > TRUE AVSK, SA = 10°. 
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BIGURE 19. FILTER OUTPUT FOR EST. FOS F(DET), 
EST. AVSK > TRUE AVSK, SA = 10°. 7th TRUE 
POINT SUPPRESSED BY 8th FILTER OUTPUT 
POINT. 
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FIGURE 20. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F(DET), EST. AVSK < TRUE AVSK, AND 
Sara 5e: 
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FIGURE 21. FILTER OUTPUT FOR EST. FO HIGHER OR LOWER 
THAN F(DET), EST. AVSK# TRUE AVSK, AND 
SA = 5°, 
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FIGURE 22. FILTER OUTPUT FOR EST. FO< F(DET), EST. 
AVSK > TRUE AVSK, SA = 5°. 
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BIGURE 23, FILTER OUTPUT FOR EST. FO> F(DET), EST. 
AVSK * TRUE AVSK, SA=5. 
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FIGURE 24. FILTER OUTPUT FOR EST. FO> F(DET), EST. 
AVSK =2-TRUE AVSK, SA = 5%. 
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~V. CONCLUSIONS 


The automatic frequency tracker developed in this report has been 
shown to be adequate for tracking a doppler-shifting frequency ina 
favorable signal-to-noise ratio environment. The width of the tracking 
oomindow.would:have.to he.increased for a higher speed target unless the 
transforms were computed at closer time intervals. A loss of signal 
would result in the last center bin frequency being retained as the signal 
being tracked unless a spurious signal with more power than the chosen 
minimum appeared in the window. Therefore, under marginal signal- 
to-noise ratio conditions, an operator would be required to monitor the 
frequency tracker output to insure that after a "lost track'' condition 
occurred, any "regain track"' operation was indeed occurring on the 
previous signal being tracked. 

The procedure of using Equation (3) to extend the frequency resolu- 
tion of the output above that of the FFT processing is dependent upon 
having a reasonable approximation of the ambient noise level. Comput- 
ing a new estimate in each transform results in higher power levels 
than the true ambient noise level being used when the target is near CPA. 
Since the ratio of two numbers, A/B, is not equal to the ratio (A-K)/ 
(B-K) where K is a constant, the subtraction of the corrupted estimate 
of ambient noise level causes an error in the value of the resolved 


frequency. However, in comparing the amplitudes of the signal and 
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noise estimates around CPA (from Tables I, II, and III), the signal is 
150 to 200 times as large as the corrupted ambient noise estimate. The 
resultant error on the resolved frequency is on the order of magnitude of 
10-4Hz and was considered negligible. As the range at CPA increases, 
the signal power level decreases and its corrupting influence on the 
ambient noise level approximation is decreased. Thus as long as the 
increased average power level in the window is caused “ the target 
being tracked, no appreciable error occurs and the resolution of the 
overall system has been extended beyond the resolution of the FFT 
processing. 

With real frequency inputs and reasonable initial estimates of the 
rest frequency and target speed, the Mitschang Filter provided target 
tracks in excellent agreement with the known emitter locations and 
probable emitter locations. Thus a true single sensor fix appears 
possible using existing hardware and techniques. Further testing is 
needed where exact target locations are known and real measured bear- 


ings are available corresponding to the real frequency measurements. 
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